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Abstract 



The Milky Way is made up of a central bar, a disk with embedded spiral 
arms, and a dark matter halo. Observational and theoretical constraints for 
the characteristic parameters of these components will be presented, with 
emphasis on the constraints from the dynamics of the Milky Way gas. In 
particular, the fraction of dark matter inside the solar radius, the location 
of the main resonances, and the evidence for multiple pattern speeds will be 
discussed. 



1 Introduction 



In our current understanding of the universe galaxies evolve within dark 
matter halos during the expansion of the universe. Very successfully, the 
theory of hierarchical clustering and merging of dark matter halos explains 
the structure on large scales and the variety of galaxy masses and types 
observed. Yet, on smaller scales the details of dissipation, star formation, 
and secular evolution, which also depend on non-gravitational forces, be- 
come more important and need to be addressed. Observations and detailed 
modeling is required to help understand the processes of galactic evolution. 

Our own galaxy is particularly important in this respect, because it 
allows much more detailed observations, down to stellar proper motions 
around the central black hole and planet formation around young stars. 
Many new large surveys of the visible matter inventory of the Milky Way are 
under way: the survey RAVE hopes to measure radial velocities and chem- 
ical abundances of the brightest 50 million stars, the space interferometry 
mission SIM will accurately determine the distance to the galactic center and 
the local standard of rest (LSR) velocity, and GLIMPSE, a Spitzer Legacy 
Science Program, already built a catalog of 47 million sources in the longi- 
tude range 10° < |/| < 64° within the galactic plane |6| < 1°, to be extended 
to fill the gap within \l\ < 10°. The goal of these su rveys is to understand 
the Milky Way's present state and for mation history ( Benjamin et al. 20051 : 
Maiewski et oilEool ISteimnet^EB . 



Detailed modeling of the Milky Way is required to make sense of these 
data sets. Because of our position within the galactic plane, we cannot 
grasp the large-scale structure easily. With certain symmetry assumptions, 
however, it is possible to recover the present 3-dimensional structure. In the 
following, we review the steps necessary for deriving such a model. 

2 Luminosity model 

Traditionally, mass models for the Milky Way have been constructed from 
parametrized analytical descriptions of the individual components: bulge, 
thin and thick dis k, bar, and h alo. The most successful axisymmetric model 
published by iKentJ <jl992h who used IRAF data to describe the stel- 
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lar matter distribution. Later bwek et al\ l|l99fj ) used the COBE/DIRBE 



near-infrared maps of the inner Galaxy to fit a tri-axial bulge and double- 
exponential disk to the near-IR data. This approach was so successful - 
for the first time the bar of the Milky Way was directly seen — that a 
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number of more detailed studies were triggered. However, any division into 
components is to some degree arbitrary, and parametric models have serious 
limitations. 

Another route was taken by Fuxl ( 1997ft . who generated a large number of 
bar forming n-body models. By selecting the model closest to the observed 
COBE/DIRBE light distribution he was able to reproduce the observations. 
The method is limited by the number of models one can generate; and there 
is no guarantee that and how well the Milky Way fits in the space spanned 
by the n-body models evolved from the chosen initial conditions. 

A more generic approa c h was adopted by Binnev and Gerhard! (|l996ft 
and Bissant z and Gerhardl ((20021 ) . who developed and applied new non- 
parametric methods, to fit the observed near-IR light distribution as mea- 
sured by COBE/DIRBE. Both methods require the specification of an initial 
model which is subsequently adjusted to simultaneously minimize the de- 
viations from the observations, smoothness and penalty term. With some 
reasonable symmetry assumptions, these methods yield the approximate 3- 
dimensional light distribution of the inner Galaxy on a model grid which 
can then be analyzed to extract parametric information, such as disk scale 
length or bar axis ratio and radial extent. 

When adding a gaseous component, the bar in this model can drive spi- 
ral arms similar to the observed spiral patter n. However, a more realistic 
treatment of the spiral arms was introduced by Bis santz and Gerhardl ()2002ft 
who allowed for deviations from triaxial symmetry, and added an analytical 
spiral arm model known to fit the observed spiral arm pattern rather well, 
and took it into account when inverting the light distribution. Interestingly, 
a model with spiral arm contribution to the COBE/DIRBE light makes the 
deduced bar longer and more elongated compared to models without spiral 
arms. Furthermore, only a model with spiral arm contribution can pro- 
duce the observed asymmetr ies in the red clump line-of-sight distributions 
( Bissantz and Gerhardl (2002ft . This result demonstrates how important it 
will be to combine different data sets from different observations. 

The modeling results depend to some extent on the assumed solar galac- 
tocentric radius R® and local standard of rest (LSR) velocity V@. iReidl 
(1993) reviewed measur e ment s of Rq and found 8 ± 0.5 kpc, while more re- 
cently E isenhauer et al. I (|2003h estimated the distance to the galactic center 
black hole to be 7.94±0.42kpc, still with a relative error of 5%. The LSR ve- 
locity is then inferred from H ipparcos data, coming out to be 218 ±21 km s _1 
( Feast and WhiteloolJ Il997h . Most of the uncertainty comes from the Rq 
error. A new analysis, t aking into account p reviously ignored systematic ef- 



fects, was published bv lKalirai et all (j2004T ) showing that the LSR velocity 
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might be lower, only 203 ± 25kms 



-i 



Reid and Brunthalerl Ho04) 
other hand found a higher value from the proper motion of Sgr A*: V© = 
236 ± 15 km s _1 . Here we assume a solar galactocentric radius Rq = 8kpc 
and a LSR velocity Vq = 220 km s -1 . 



3 Mass model 

Assuming a spatially constant mass-to-light ratio, the NIR luminosity model 
can be converted to a mass model. In reality, the mass-to-light ratio of a 
stellar population depends on age and metallicity of the underlying stellar 
population, which would vary within the Milky Way, but these dependencies 
are weaker than in the optical. The luminous Galaxy is probably embed- 
ded in a dark matter halo, which makes the circular rotation curve flat in 
the outer parts. Combining the luminosity model and the constant M/L 
assumption for the stars with an analytical model for the dark matter dis- 
tribution, a mass model can be constructed with only a few free parameters: 
the mass-to-light ratio for the visible component(s), the position angle of 
the bar, the pattern speeds of bar and spiral pattern, and the halo model 
parameters. 

For the invisible dark matter halo, we assume a logarithmic density pro- 
file with two parameters: scale radius a and velocity V. We find below that 
the details of the assumed profile do not matter, as long as the rotation 
curve of the dark halo is proportional to radius within the region probed 
by the data, i.e. inside Rq. Therefore, the dark matter halo model truly 
has only one free parameter, the slope at small radii V/a, where a can be 
anything larger than lOkpc. 

The combined total rotation curve can be compared to the observed 
terminal velocity curve of the gas. When we assume that the luminous 
matter in the Milky Way has the maximal M/L allowed by the terminal 
curve, and a dark matter halo scale radius of a = 10.7 kpc, the best fit 
is obtained with V = 235 km s -1 . The result can be seen in Fig. Q The 
deviation for r < 2.5 kpc is due to the bar which forces the gas on non- 
circular orbits. At larger radii, the agreement is very good, and apart from 
two small humps at 4 kpc and 6 kpc, the model fits as well as the scatter 
in the data allows. When the orbital motion of the Sun is added to the 
model's terminal velocity curve, we obtain the rotation curve in Fig. |2I One 
sees that the contribution of dark matter inside Rq in the model is less than 
the contribution of visible matter. 

From Fig ^ we see that a surprisingly good fit to the data is obtained 
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Figure 1: Comparison of model (solid curve ) to a collectio n of terminal 
veloc ity curve data (symbols: 1 2 CO data f rom Clemens! (1985), and HI data 
from Burton and Lisztl ( 19931 ): iFich et al\ dl98fll )). The lower dotted curve 
corresponds to the visible light contribution. For this fit, we assume Rq = 
8 kpc and solve for the LSR velocity. The fit is required to match zero 
terminal velocity at Rq and the only free parameters varied are the halo 
velocity scale V and the mass-to- light ratio for the visible matter, which 
depends weekly on the assumed V. By adjusting V we were able to match 
the observed terminal velocity curve; the LSR velocity for this fit is Vq = 
221 km s -1 , in agreement with observations. 
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Figure 2: Rotation curve for dark matter (dotted), visible matter (dashed), 
and total (solid) for the model shown in Fig. ^ 
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when the contribution of visible matter is maximized and a constant mass- 
to-light ratio is assumed. Additional evidence for this result comes from 
microlensing statistics obtained for selected fields within the bulge. The 
luminous mass model is only m arginally able to reach the high microlens- 
ing o ptical depth measurements ( Bissantz et al.lll997 ; Bissantz and Gerhard! 
2002). Increasing the dark matter component in the model would lower the 



model's predictions of the microlensing optical depth, and thus increase the 
discrepancies with the experimental results even more. 

Turning the analysis around, we can assume that we got the contribu- 
tion of the visible matter to the rotation curve right, and then calculate the 
rotation curve of the dark matter halo, i.e. what is left for the dark matter 
after subtraction of the known visible matter (see Fig. . The data up to 
Fin ar e taken from Fig. and those beyond Rq are taken from Merrifieldl 
( 19921) 1 who derived them from the HI layer thickness. Clearly, within the 
solar orbit, the dark matter halo contribution to the rotation curve rises ap- 
proximately linearly with radius. At small radii, the data again are strongly 
affected by non-circular motions in the bar potential. The data seem to 
suggest that we start to see the turnover, and that the halo scale radius 
is ~ 15 ± 5kpc, but the scatter in the dat a points beyond Rq is quite 
large. Moreover, iBinnev and criticized that the increased 

scatter and apparently linear rising rotation curve beyond the solar radius 
can be explai ned naturally if trac ers beyond Rq are concentrated in a ring 
at 1.6-Rq. If Binnev and Pehnenl are correct, this would imply that lower 
values for the halo scale radius a are also possible. 

The analysis so far has assumed that the gas motions can be approxi- 
mated by circular orbits. However, in the bar region (R < 3kpc) this is not 
the case and at larger radii we may have significant perturbations due to 
spiral arms or orbit family changes at resonances. A more detailed analysis, 
however, is possible with gas dynamical models which take these factors into 
account. 



4 Orbits and resonances 

The most important orbital families in the galactic plane of a barred galaxy 
are the so-called x\- and X2-orbits. The former exist inside the corotation 
resonance (CR), where one rotation of the bar takes as long as one circular 
orbit period. The £i-orbits are shaped by the periodic forcing of the bar and 
are in phase and elongated with the bar. The X2-orbits are perpendicular to 

1 after scaling to the values of Rq and Vq assumed here. 
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Figure 3: Dark matter halo circular rotation curve deduced from 
the terminal velocity curve after subtraction of the visible mat- 
ter contribution (all symbols). The series of linearly rising solid 
curves represent isothermal dark matter halo models with parameters 
(a,V) € {(10,223), (15,316), (20,412), (25, 509)}(kpc, kms" 1 ) and V/a ~ 
21 kms -1 kpc -1 . All these models fi t the data inside Rq = 8 kpc equally 
well. Beyond Rq, we show data from iMerrifieldl (^992) scaled for the values 
of Rq and V© used here. This figure demonstrates that data beyond Rq is 
crucial for distinguishing between different halo mass models. The rotation 
curve for the visible matter is overlaid for comparison (dotted line). 
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the bar and typically almost round. They exist between the inner Lindblad 
resonance (ILR) and the center or between the two ILRs if present. Lind- 
blad resonances happen when the epicycle frequency is twice the difference 
between the pattern speed and orbital angular speed. Inside corotation we 
can have typically one or two ILRs; outside corotation there is always a 
single outer Lindblad Resonance (OLR). Between CR and OLR, orbits are 
again elongated perpendicular to the bar but almost circular in shape. 

The xi-orbits are considered main building blocks of the stellar bar, 
because they support the bar potential. They are als o responsible for th e 
sharp peak in the Milky Way's terminal velocity curve ([Binnev et aZ.M l991). 

Using the atomic hydrogen and molecular gas emission all sky surveys, 
which not only provide terminal velocity data for the inner galaxy, but also 
information about spiral arm tangents, etc., we can use hydrodynamical 
models to find further constraints for the mass model. 



5 Modeling the gas dynamics 

Given the mass model, one can find quasi steady state gas flow solutions 
by assuming some equation of state for the gas and solving the Euler equa- 
tions. Here, the equation of state is assumed to be isothermal with a sound 
speed of lOkms -1 ; and for the modeling we use a 2-dimensional smooth 
particle hydrodynamics (SPH) code. Artificial viscosity in the code is used 
to simulate shocks properly. The background potential is split into three 
components: axisymmetric, bar, and spiral arms. The bar and spiral arm 
components rotate with fixed pattern speeds Sl^ and Sl s . Initially, the gas is 
set up as an exponential disk on circular orbits in the axisymmetric part of 
the potential. The non-axisymmetric forces due to bar and spiral arms are 
introduced slowly in the model to avoid transient features. The resulting gas 
flow is stationary in the rotating frame of the bar if bar and spiral pattern 
speeds are identical, or periodic if they differ. For a given set of model pa- 
rameters, longitude-radial velocity diagrams for the gas are then computed 
as would be measured by an observer rotating with the LSR, from which we 
obtain the model terminal velocity curves. 

For the case of a single pattern speed, the model is evolved until the 
gas flow is almost stationary 2 in the bar frame of reference. For the case 
of separate pattern speeds for the bar and spiral patterns, the time step 
corresponding to the gas flow closest to the observed gas flow is found 

2 A small evolution with time is always present because of the mass inflow towards the 
center. 
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by comparison with th e obse rved atomic and molecular gas kinematics. 
Englmaier and Gerhard! (|l999l ) have shown that the terminal velocity curve 
can then be well understood in the context of gas on orbits typically found in 
barred galaxies. Their gas dynamical model in the barred Milky Way model 
potential fits the observed terminal velocity very well (see Fig. QJ). Note, 
that the high peak at I = 2° is explained by non-circular motions in the bar 
potential and is matched better in models with higher numerical resolution. 
The whole process is iterative, because changing the position angle of the 
Sun with respect to the bar major axis requires recomputation of the mass 
model with the corresponding symmetry assumption, and recomputation of 
the gas flow in the modified potential. 



6 The galactic bar and the spiral arms 

The best near IR model of Bissantz and Gerhard! |2002h has a 3-dimensional 
bar with axis ratios 1 : b : c where b = c = 0.35 ± 0.05. The angle between 
the bar major axis and the line from the Sun to the Galactic Center is 20°, 
with the bar's near end at positive longitudes. 

The corresponding mass distribution, however, contains no direct in- 
formation about dynamical parameters such as pattern speeds. If we can 
assume that the Milky Way's bar and spiral patterns each rotate with a 
constant angular speed for at least a few rotations, and that the gas flow 
reaches a quasi-steady state, then a comparison of the model with the ob- 
served gas flow can yield constraints on the current values of both pattern 
speeds in the Milky Way. 

In the gas dynamical model we observe spiral arms that cause large 
deviations from circular velocity within corotation and only small deviations 
outside corotation of the bar. In the Milky Way, there is the so-called '3- 
kpc-arm' which extends to 3kpc in radius and has strongly non-circular 
motions. On the other hand, the so-called molecular ring material between 
4kpc < r < 7kpc appears to be on nearly circular orbits. Thus we can 
place corotation between those limits; in the following we set the corotation 
radius to 3.5 kpc. 

From observations and models of other barred galaxies it is known that 
the corotation radiu s is typically 1.2 ± 0.2 times the radius of the bar ends 
( Athanassoulalll992h. Larger values would b e possible if the bar were able to 
slow down (jPebattista and Sellwoodll200of ). Our corotation derived above 
is consistent with the condition i?cR/-Rbar = 1-2 ± 0.2, however, if the bar 
length turns out to be underestimated, we would need to increase the coro- 
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Figure 4: Gas dynamical model results for the Galactic termi- 
nal velocity curve. Data p oints are from various observations (see 
Englmaier and Gerhard! (|l999h for details). The upper (lower) part of the 
diagram shows data and models for positive (negative) galactic longitudes. 
The solid curve is the best model including dark matter. The dash-dotted 
curve is a model without dark halo. 
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Figure 5: Gas response for models with spiral arm pattern speed set to 60 
(left), and 20kms _1 kpc~ 1 (right). 



tation radius accordingly. Dehnenl |2000) found evidence from local stellar 
st reams that the OLR of the bar is located very close to the Sun, according 
to Muhlbauer and Dehnenl (J2003) . at ~ 0.92i? Q . If so, we can calculate the 
corotation radius for our model and find Rcr = 4.0 kpc. This is in good 
agreement with our previous estimate. 

The gas dynamical model in the barred mass model without spiral arms 
for ms a 4-armed spiral patte r n sim ilar to the one deduced from H H-re gions 
by Georgelin and Georgelin dl976h: see lEnglmaier and Gerhardl ( 1999h for 
the comparis on. IBissantz and Gerhard! (|2002l ) used the more realistic spiral 
arm model of O ritz and Leoin ej |l993F ) to improve the inversion of the inner 
Galaxy, and Bissantz et al\ (|2003 ) found the approximately periodic gas flow 
in this modified potential, assuming separate pattern speeds, for the bar and 
spiral arms. 



7 Multiple Pattern Speeds 

When we allow two pattern speeds in the model, i.e. the spiral pattern 
is given a different pattern speed from that of the bar, then the gas flow is 
changing in a periodic fashion with the beat frequency £1^ — O s . This periodic 
oscillation is strongest in the transition region between bar and spiral pattern 
around the bar's corotation. Because of the lower pattern speed, the spiral 
pattern extends to larger radii; its outer radius is approximately given by 
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Figure 6: {l,v) diagrams for models with fast (60kms _1 kpc _ ; top) and 
slow spiral arms (20kms -1 kpc -1 ; bottom) as shown in Fig. El Only the 
slow model shows low density voids (marked as 'V') similar to those known 
from observations. 
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its OLR. 

Dist ant faint gas eous spiral arms have b een found out to radii 3 of 12 - 
16 kpc dDavieslll972h . and even 17-22.6 kpc (|McClure-Grimths et a,lhm& \. 
However, it could well be possible that the outermost spiral arms are tran- 
sient and do not follow the same pattern speed. Hence, it is not practical to 
use the radial extent to estimate the spiral arm pattern speed. 

In the model with two pattern speeds, most of the time the spiral pattern 
beyond bar corotation smoothly connects to the inner spiral pattern driven 
by the bar (Fig. 03). However, because of the different pattern speeds, the 
outer spiral arms regularly separate from the bar-driven inner arms around 
corotation. On the other hand, when spiral pattern and bar rotate with 
the same pattern speed, the spiral pattern weakens in the corotation region 
and disappears. This is consistent with density wave theory, which asserts 
that density waves become evanescent in the region around corotation, and 
decay. When two pattern speeds are used, we observe in the model that the 
spiral pattern passes through the bar's corotation and through the spiral 
arms' corotation radius without weakening. A detailed comparison of the 
model with the observed gas flow shows that certain gas poor regions in the 
observed longitude-radial velocity diagram correspond to inter-arm regions 
in the gas flow near CR which only exist in the corotation region when spiral 
arms and bar rotate with different pattern speeds (Fig.^J). The observed gas 
dynamics clearly shows such gas poor regions, hence the spiral arms must 
rotate more slowly than the bar fBissantz et a/.l Eo03). Approximate agree- 
ment with observations is found when a pattern speed of 20kms~ 1 kpc -1 is 
assumed for the spiral arms. 



Detailed observations by iNaoz and Shavivl (|2f)05fi of the spatial separa- 
tion of spiral arms and young star clusters that have moved away from their 
spiral arms, also indicate not only that the spiral pattern rotates much slower 
than the bar, but also that multiple pattern speeds are present in the spiral 
pattern which fall into two groups around 18 and 30kms _1 kpc -1 . This 
observational result is also in agreement with the theory of Lin & Shu (1964) 
and a general result from n-body simulations of disk galaxies, namely that 
spiral arms are in fact not steady -state but recurrent and multiple modes 



overlap each other (e.g. lF^ll99 



3 Rescaled to correct for the assumed distance of 8 kpc to the Galactic Center. 
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8 Discussion 



In this review, we gave an overview of the results obtained from reconstruc- 
tion of the present state of the Milky Way using near-IR luminosity and 
gas kinematics complimentary data sets. From the luminosity data alone, 
one could not infer anything about the dark matter component. Likewise, 
the observed gas kinematics alone is difficult to interpret and has lead his- 
torically to wrong conclusions, such as steep drops in the inferred rotation 
curve, large dark matter components, and explosions in the galactic center. 
We avoid using parametric models for describing the important bulge region 
which is best described as a barred bulge. This bar in the galactic center 
explains the observed asymmetry in the near-IR light distribution within 
the bulge region as well as the high non-circular gas motion inside 3 kpc 
galactocentric radius. 

The inferred mass density is in good agreement with micro-lensing op- 
tical depth observations which provide an important independent test. The 
molecular ring material, is interpreted as four tightly wound spiral arms 
in approximate circular rotation outside the corotation radius of the bar. 
Models indicate that the spiral arm pattern rotates much slower than the 
bar, allowing the arms to extend beyond the solar orbit. Most imporantly, 
the model has provided evidence for the maximum disk hypothesis and the 
amount of dark matter present inside the solar radius. Future and presently 
ongoing observational surveys will provide even more detailed data which 
will allow us to improve the model further. The next revision of the Milky 
Way models will probably allow more detailed studies of the spiral arms and 
locate more resonances. 
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